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Intrinsic heavy quarks in hadrons emerge from the non-perturbative structure of a hadron bound state [l| and are a 
rigorous prediction of QCD [1,Q. Lattice QCD calculations indicate significant intrinsic charm and strangeness prob¬ 
abilities 0, Q . Since the light-front momentum distribution of the Fock states is maximal at equal rapidity, intrinsic 
heavy quarks carry significant fractions of the hadron momentum. The presence of Fock states with intrinsic strange, 
charm, or bottom quarks in hadrons lead to an array of novel physics phenomena il- Accurate determinations of 
the heavy-quark distributions in the proton are needed to interpret Tevatron and LHC measurements as probes of 
physics beyond the Standard Model [10. Determinations iEll [T^ of the momentum fraction carried by intrinsic 
charm quarks in the proton typically limit (a:)ic ~ 0(1%) at 90% CL, consistent with the EMC analysis of their 
charm structure function data |13l| and the large rate for high-p^ pp c'yX reactions at the Tevatron 141; however, 
a precise determination of (a:)ic has proven elusive. 

The recent letter by P. Jimenez-Delgado et al. (JDHLM)J0 reports (x)ic = (0.15 ± 0.09)%. The authors include 
low-energy data from the ed (p) —>• e'X SLAC experiment [la| in their global fit. They find that this data set places 
strong constraints on intrinsic charm, although, by their count, only 157 of 1021 data points have in excess of the 
charm hadronic threshold: ~ 16GeV^. The SLAC measurements of ed{p) e'X have an overall normalization 

(systematic) error of ± 1.7 (2.1)%, and a relative normalization error of typically ±1.1% 16|. The SLAC data points 
in the germane regime of LF^ > 16 GeV^, x > 0.1 have even larger statistical uncertainties. 

It is clearly challenging to identify the contribution from charm quarks to the inclusive structure function if only the 
scattered electron is detected. In addition to the valence and sea-quark distributions, there are other contributions 
to the inclusive cross section which need to be determined to high accuracy in order to discern the intrinsic charm 
component at the level claimed by JDHLM; this includes higher-twist corrections at high x, the strange and bottom 
quark contributions, as well as the accurate implementation of the suppression of charm production at threshold and 
nuclear target effects. Their analysis depends on uncertain parameters and theory assumptions. For example, JDHLM 
model higher-twist effects as an isospin-independent, phenomenological multiplicative factor on top of a leading-twist 
structure function with target-mass corrections 0. Their higher-twist model does not include enhancements at x ^ 1 
resulting from hard subprocesses such as e[qq\ e'qq where multiple quarks interact 18, 0|. Such processes depend 
strongly on the diquark charges, making them quark flavor (isospin) sensitive, and they contribute in the same large 
X domain as the charm signal. In addition, the target mass-corrected structure functions used by JDHLM remain 
nonzero as x ^ 1 0 which is problematic, and they neglect intrinsic strange and beauty quark contributions 0,0. 

JDHLM assess their PDF errors using a tolerance criteria of = 1 at Icr; however, the actual value of Ay^ 
depends on the number of parameters to be simultaneously determined in the fit — their assessment of a single 
parameter error requires that the other parameters be fixed at their values at the global y^ minimum 211. JDHLM 
report (x)ic = (0.15 ± 0.09)% 0 corresponding to Ay^ = 1 at Icr (68% CL) and also (x)ic ~ 0.5% at 4cr. In 
order to set a 4cr limit, all of the other parameters must be varied so as to yield a minimum y^ as the parameter of 


interest is changed [22l - l25|| . We note Refs. [26j, [ 23 , e.g., contain 25 PDF parameters in leading twist, and Ref. [27 1 
contains 12 more higher twist parameters. Since one would expect nontrivial correlations and near degeneracies in 
a many-parameter fit, the apparent agreement of the 4 x Icr assessment with the 4cr limit suggests that the other 
parameters were not properly varied as (x)ic was changed, making the reported limit too severe. 

It is clear that the SLAC single-arm measurements cannot unambiguously identify an intrinsic charm contribution 
to the nucleon structure function even at the 1% level because of statistical and systematic uncertainties, both 
experimental and theoretical. JDHLM have excluded the EMC measurements of the charm structure function [0 . 
citing a “goodness of fit” criterion. However, statistical criteria alone cannot exclude data sets. The fit to the EMC 
iron target data would be improved by including the QCD nuclear effects described in Refs. 28. 0. 

In summary, JDHLM claim that the momentum fraction carried by intrinsic charm quarks in the proton is (x)ic = 
(0.15 ± 0.09)%; we do not find this conclusion warranted. 
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